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ABSTRACT 



A method to construct Hamiltonian theories for systems of both ordi- 
nary and partial differential equations is presented. The knowledge of a 
Lagrangian is not at all necessary to achieve the result. The only ingredients 
required for the construction are one solution of the symmetry (perturbation) 
equation and one constant of the motion of the original system. It turns out 
that the Poisson bracket structure for the dynamical variables is far from be- 
ing uniquely determined by the differential equations of motion. Examples 
in classical mechanics as well as in field theory are presented. 
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Hamiltonian methods are widely used in connection with problems in 
classical, quantum and statistical mechanics, fluid dynamics, optics, solid 
state, molecular, atomic, nuclear, particle and plasma physics, in both clas- 
sical and quantum field theoretical systems. Quantization schemes as well 
as group theoretical symmetry methods are examples of subject matters in 
which Hamiltonian structures are useful. Hamiltonian theories are usually 
constructed starting from the knowledge of a Lagrangian, by well established 
methods for both the cases of regular and singular Lagrangians [|l], ^, ^, 
For different reasons, one may try to quantize or to construct Hamiltonian 
structures for classical systems of differential equations, without recourse to a 
Lagrangian [@, |^, . Several authors have been successful in creating Hamil- 
tonian theories from scratch for different examples, mostly in fiuid dynamics 
(an excellent review is presented in [Q]) and in field theory |p, but no general 
method seems to exist for constructing a Hamiltonian structure starting from 
the equations of motion only, without using at all either the explicit form or 
the existence of a Lagrangian formulation for the system at hand. 
The purpose of this note is to present a general technique to construct Hamil- 
tonian theories starting from the equations of motion, using one symmetry 
transformation and one constant of the motion, without recourse to the La- 
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grangian of the system of equations, which may even fail to exist. A com- 
pletely different approach has been used to construct a quantum model for 
a non-Lagrangian cosmological model in ITOf); and infinitely many Hamil- 



tonian structures of the spinning top in |Tl|. Let us first define what it is 

usually meant by a Hamiltonian theory. 

Consider an autonomous first order differential system, 

fa/ 6 



/V) a,b=l, ,N. (1) 



dt 

Of course, a differential system of any order can be easily cast in first order 
form by defining extra variables in the standard textbook fashion. A Hamil- 
tonian structure for ([1|) is defined in terms of an antisymmetric symplectic 
matrix J"-^{x'^) and a Hamiltonian H{x^) which satisfy 

jab^_jba ^^^^^^ ^2) 

r\d J'^ + J^\d J'^ + J'^^ = , (3) 



and. 



r)T-f 

r^T^^'f". (4) 



The Poisson bracket between any two dynamical variables A[x°') and B{x^) 
is defined by 
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and it satisfies all the usual algebraic and differential properties as it can be 
easily inferred from the antisymmetry condition (0), the Jacobi identity (|^), 
and the definition (H). 
If, in addition, 

det 7^ , (6) 

is required, then the symplectic matrix is regular. It is important to remark 
that, sometimes, condition (|^) cannot be met. If N is odd, as it happens 
for Euler's equations, (H) is never satisfied because of (0) as it can be seen. 



for instance in |TT| and in one of the examples below. In Dirac theory, 
condition (|^) for Dirac brackets is not satisfied because there are dynamical 
entities, called second class constraints, which have vanishing Dirac brackets 
with any variable 0]. In fluid dynamics, the functions which have vanishing 
Poisson brackets with any dynamical variable are called "Casimir functions" , 
inspired on the well known group theoretical terminology for operators which 
commute with any element of the group |12[. Therefore, it is convenient 



to adopt a flexible attitude regarding condition (||), and take conditions (Q), 
(^, and (Q) as defining a Hamiltonian theory. Of course, the usual textbook 
Hamiltonian structures satisfy all of them. 

It is straightforward to prove that conditions @j and imply that 



is a time independent constant of the motion of system (0) defined by the 
condition 

^^r-C,H^,, (7) 

which can be equivalently stated by saying that the Lie derivative of H along 



/ vanishes. A brief account on Lie derivatives may be found in |]T3 . 

Let us now derive the symmetry (perturbation) equation of (|l]). (For a 

detailed discussion, see [|14|)- Consider the transformation 



F = x° + er]\x\t) , (8) 

where e r]°'{x'^,t) is a small perturbation which maps solutions of (|l|) in so- 
lutions of the same equation, up to first order in e. The equation that the 
perturbation vector rj satisfies is, 

dt rf + r)\, f - = , (9) 

or, 

(dt + Cf)v^ = . (10) 

It is not difficult to prove that K, the deformation of H along i], defined by 

r) J-f 
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is also a constant of the motion for the same system, if 77 satisfies Eq. (|T0|). 
By the same token, a new symmetry transformation f] which satisfies Eq. 
(p!oD can be constructed using a symmetry transformation rj and a constant 
of motion K by 

= I • (12) 



A detailed account of these results may be found in [|I5 . 
Let us now compute the Lie derivative of J"-^ along / 

Cf r' = r\c r - f,c - J''" r,c ■ (i3) 

It is an straightforward exercise to prove that 

Cf = . (14) 

using Eqs. (0), (|^), and (^). Note that Eq. is not needed in the proof. 
Therefore, a symplectic matrix must have vanishing Lie derivative along /. 
Nevertheless, this condition is not sufficient to fullfill simultaneously the re- 
quirements (0), (0), and (^. To construct a symplectic matrix J""^, let us start 
by considering an antisymmetric matrix according to the following ansatz 

jab = far,' - f rf , (15) 

where 77 satisfies (H) and has been normalized using (|T^) in such a way that 
J"'' fulfill (^ identically. Of course, condition (0) is trivially met. The Jacobi 
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identity (Q) imposes the following condition 

Cf T]" + J™ Cf rf + Cfr]^ = 0, (16) 

which is satisfied by a particular, time independent symmetry vector rjo which 
solves (p!OD defined by 

dt iio" = -CfTio'' = . (17) 

A more interesting solution r]i is given by the condition 

dt vi' = -Cfv^'' = Xr , (18) 

which will be most useful in many instances. Note that both solutions pro- 
duce symplectic matrices with vanishing Lie derivatives along /. 
We have have thus contructed a Hamiltonian structure for (|I|) based on the 
knowledge of just one symmetry vector (either rjQ or r]i) and only one con- 
stant of the motion, H, of the system under consideration (assuming a non 
vanishing K, which can be easily achieved as it will be seen in the examples). 
A few comments seem in order. First, a solution similar to the one described 
by (p!8D, with the Lie derivative of the symmetry vector along / proportional 
to the symmetry vector itself, although it satisfies the Jacobi identity, is in- 
compatible with (^). Second, the rank of the symplectic matrix just derived 

8 



is two. Therefore, it will be, in most cases, singular. A procedure to en- 
large its rank will be described below. Third, it is obvious that the method 
we have presented will, in general, yield a Hamiltonian structure written in 
terms of non-canonical coordinates. Nevertheless, that is the most we can 
hope for in the case of non-Lagrangian systems (which are always described 
by a non-commutative geometry). This problem is dealt with in some detail 
in 1^. Fourth, even though this procedure differs from the usual one for the 
case of Lagrangian systems, it may sometimes reproduce the well known re- 
sults in terms of canonical coordinates, as it is shown in one of the examples 
below. 

Let us now consider some examples. The systems may be completely de- 
scribed by the evolution vector /, or the equations of motion written in its 
first order version. The Hamiltonian structure may be completely determined 
by one constant of motion, the Hamiltonian H, and one symmetry vector rjQ 
or r]i. Sometimes, we will need to make use of the the normalization given 
in (0). 

Example 1.— One dimensional monomial force. 
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This example is defined by the equations of motion 

f ^ x\ f ^ -c{n + l){x'r, (19) 
while a Hamiltonian 

^ = ^ + c (x^)"+^ , (20) 
and one symmetry transformation are given by 

- + tf\rf^ x^ + t f . (21) 

This is, of course, a very trivial example, which, nonetheless, shows how the 
scheme presented here can reproduce the usual results. In this case, the sym- 
plectic matrix is regular. The harmonic oscillator and the free particle are 
special cases in this example. Note that this treatment can be extended to 
any number of dimensions provided the force be a homogeneous function of 
degree n in the coordinates. 

Example 2.- Euler's Top. 
Consider the equations of motion of Euler's top 

— ^-e^^%L^=f i,j,k^l,2,3 , (22) 
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with 

L 



^j = -r > (23) 

where Lj = L* and Qi = fi* are the components of the angular momentum 
vector and the angular velocity vector in the i*'^ principal direction respec- 
tively and the Jj is the eigenvalue of the tensor of inertia of an asymmetrical 
top along the i**^ principal usual. 

Let us now look for symmetries of the equations of motion. With this purpose 
in mind, multiply the angular momentum by some constant factor A. This 
operation introduces a factor in the right hand side of the equation of mo- 
tion (|2^). The same result is achieved in the left hand side of the equation 



if, in addition, time is multiplied by the inverse factor A^^. These opera- 
tions performed simultaneously constitute a finite symmetry transformation 
for (|2^) . One can deal with an infinitesimal version of it by considering 
A = 1 + C infinitesimally close to one, to get the transformation 

6U = QL\ 6t = -Ct , (24) 

Note that a transformation such as (|2^) may equivalently be written as (see, 
for instance [|14|) 

rj^ = aL' + te'^''n,Lk) , (25) 
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leaving time invariant. It is now a straightforward matter to check that the 
transformation defined by (p5D is, in fact, a symmetry transformation for (^) 
because it satisfies @ . Note that rj^ also satisfies (|T8|) , and consequently may, 
in principle, be used to define a symplectic matrix according to (|T5|). 
It is well known that Ci and C2 given by 

C, = {L'f + {L'f + {L'f , (26) 

and 

are constants of the motion for the dynamics generated by (0). We have 
already seen that the Hamiltonian for any system must be a constant of the 
motion. Therefore, Ci and C2 are, in principle, possible Hamiltonians for the 
top. 

The deformations of Ci and C2 along rj^ do not vanish, in fact, 

^1 - ^ = 2C^i > (28) 

and 

K2 ^ ^ V' = . (29) 

We have thus found two inequivalent Hamiltonian formulations for the top, 
defined by the symplectic matrices Ji-' and J2*'' and the Hamiltonians Hi 
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and H2 given by 

Ji'' -^^UW - f ) , (30) 

Hi ^C, , (31) 

J2'' {frf - f ) , (32) 

and 

H2 . (33) 

Note that the choice of a Hamiltonian H as an arbitrary function H — 
H{Ci,C2) is also possible provided the proper normalization factor K is 
used in the symplectic matrix J'-' . In this way, we have constructed infinitely 
many Hamiltonian structures for Euler's top. 

Example 3.- Radial forces. 
This example considers non-potential radial forces defined by 

f = x'+^ , = F x\ F = F{f^,r ,r-'r), i = 1,2,3. (34) 

To construct a Hamiltonian structure, the Hamiltonian may be chosen to be 
the third component of the (conserved) angular momentum vector, while the 
symmetry transformation is, for instance, a rotation around the first axis. 
This example clearly shows the ambiguity which exists to choose both the 
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Hamitonian and the symplectic matrix. 



Example 4.— Korteweg-de Vries equation. 
The equation of motion is 

ut = - u - u^xx = f . (35) 

We are now going to construct a symmetry transformation for it. Take 



any solution of Eq. (|35D and define a new set of variables u', x', and t' by 
multiplying the old variables u, x, t by factors A~^, A, respectively. This 
operation simply produces an overall A^ factor in the equation, which means 
that the new set of variables solves the same equation which the old variables 
satisfy. We have thus constructed a finite symmetry transformation for the 
Korteweg-de Vries equation. The infinitesimal symmetry associated to it 
may be written taking A = 1 + C, infinitesimally close to one 

6u = —2Qu , 5x = C,x , 5t = 3(t , (36) 

or, equivalently, 

1] = C(-2m - X + 3t{uu^ + u^^^ ) ) , (37) 



leaving x and t unchanged |TJ]. It is now a straightforward matter to prove 
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that 1] is an infinitesimal symmetry transformation for (^) because it sat- 
isfies tlie symmetry (perturbation) equation @. We remark tliat rj satisfies 
condition (^) and therefore, it can be used to construct a symplectic matrix. 
One possible choice for the Hamiltonian density is u^. The Hamiltonian H 

H = J u^dx , (38) 

is a constant of motion, i.e., its time derivative vanishes when the usual 
assumptions about the behaviour of the fields at spatial infinity are adopted. 
In fact , the time derivative of the Hamiltonian can be written as the integral 
of a total spatial divergence (a partial derivative with respect to x, in our 
case) when the equation of motion ( p5|) is taken into account. 
One gets that the deformation K of H along r] is non-vanishing. We can 
easily see that the functional derivative of H in the i] direction is 

K = [ r]{x) dx = -3C H , (39) 

J ou{x) 

as it can be obtained by direct computation, or by taking advantage of the 
of new variables defined by multiplying the old ones by powers of A as we 
have already done above. 

Therefore, one symplectic structure is given by 

Ax,y) = -i ( fix) viy) - f{y) vix) ) . (40) 
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The equation of motion can now be written in Hamiltonian form 



ut ^[u , H] 



(41) 



where the field theoretical Poisson bracket is defined, as usual, in terms of 
functional derivatives by 



As far as we know this is a new Hamiltonian structure for the KdV equation. 
Note that other Hamiltonian densities H' can be used as well, in conjunction 
with the same symmetry vector rj, provided the corresponding deformations 
K' be used in the definition of the new symplectic matrix J'{x, y). 

Example 5.- Non-linear Schrodinger equations. 
The equations of motion are 



and its complex conjugate. One possible non-standard Hamiltonian density 
is -0*, and the symmetry vectors are 




(42) 



(43) 



V ^ - -0 - a; + 2 i ( + V'^ "0* ) 



(44) 
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and its complex conjugate. This Hamiltonian structure also appears to be 



new. 



Of course, some of these structures, have singular symplectic matrices. 
One way to increase the rank of the symplectic matrix, without altering any 
other of its properties, is the following. Assume we can find two new time 
independent symmetry vectors 772 and rj^ such that the Lie derivatives of the 
Hamiltonian H along them vanish, i.e., 

and that the Lie derivatives of r]2 along 773 as well as those of 772 and 773 along 
r]i (or 7]o) vanish. Then, the new symplectic matrix Ji"^ defined by 

satisfies all of the requirements which define a symplectic matrix (^, (^), 
and even the non-linear Jacobi identity (j^), and its rank is equal to four. 
This procedure can be repeated at will, producing an increase of two units 
in the rank of the symplectic matrix each time that it is performed. (Note 
that this construction clearly shows that the Poisson bracket structure is not 
uniquely determined by the dynamics). If, eventually, one gets a regular 
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symplectic matrix, the method presented here may constitute an alternative 
to construct a Lagrangian description of the system (|ID, yielding a novel, 
symmetry based, approach to the classical Inverse Problem of the Calculus 



of Variations [III llTl HB, [T9f 



Note that the choice of the symmetry vector (t/q or r/i) needed to define the 
symplectic matrix is determined solely by the requirement of getting a non- 
vanishing given H . We have used in the examples both time-dependent 
and time-independent symmetry vectors, to illustrate different possibilities. 
Sometimes, there may be several adequate choices of symmetry vectors for a 
given producing different Hamiltonian formulations for the same system. 
We remark that singular symplectic matrices are present in Dirac's construc- 
tion of Hamiltonian structures, as we have already mentioned above. We 
are currently investigating the possibility of applying this method, which 
naturally leads to singular symplectic matrices, to deal with gauge and con- 
strained systems, as an alternative to Dirac's method when no Lagrangian 
is available. We are also studying whether it is possible to obtain constants 
of the motion of the system at hand as Casimir functions of the singular 
symplectic matrix constructed here. 
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